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Abstract 

We have compared exact numerical results for the Lipkin model at finite 
temperature with Hartree-Fock theory and with the results of including in 
addition the ring diagrams. In the simplest version of the Lipkin model the 
Hartree-Fock approach shows a "phase transition" which is absent in the 
exact results. For more realistic cases, Hartree-Fock provides a very good 
approximation and a modest improvement is obtained by adding the ring 
diagrams. 



PACS number: 2L60.-n 



1 Introduction 



The study of hot nuclei and hot nuclear matter is of importance in heavy ion 
collisions and in supernova explosions. Theoretical treatments require the 
use of finite-temperature many-body theory which is inherently more diffi- 
cult than the corresponding zero-temperature formalism. Since in practical 
situations the many-body theory cannot be solved exactly, approximations 
are needed and it is often difficult to know how accurate these might be. It 
is therefore useful to study a model which has some realistic features, but 
which is sufficiently simple to permit an exact solution, so that the accuracy 
of various approximations can be assessed. With this objective in mind we 
shall study here the Lipkin model |l| which has been widely used at zero 
temperature. Rather little work has been carried out at finite temperature 
with the Lipkin model [Q, ^, Q and then only for the simplest version of the 
model. The principle interest in these studies was excited states and boson 
expansions. Particularly striking was the demonstration that a phase tran- 
sition can occur in the Hartree-Fock (HF) solutions. We shall point out that 
the exact internal energy shows a qualitatively different behavior. 

In addition to HF we will consider the correlations of the particle-particle 
hole-hole ring diagrams, i.e., the random phase approximation (RPA), for 
which an expression for the grand potential has recently been given ^j. 
We choose these approximations since Yang, Heyer and Kuo [0 found that 
at zero temperature HF gave a very good approximation to the exact ground 
state energy and the ring series gave a further small improvement. The 
final result was therefore very close to the true energy of the system. The 
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present study can be regarded as the continuation of the work of ref. [0 
to finite temperature. An additional motivation for studying the long range 
correlations of the ring series is that in calculations with realistic interactions 
it significantly improves the location of the saturation point in nuclear matter 
[§]. We note in passing that we shall not consider the particle- hole ring series 
since it gave a negligible effect in the zero temperature work of ref. 0. 

The organization of this paper is as follows. In subsec. 2.1 we discuss our 
method of obtaining exact numerical solutions for a slightly modified version 
of the original Glick, Lipkin and Meshkov model [|l| . The approximate many- 
body methods with which we compare, namely HF and RPA, are discussed 
in subsec. 2.2. The comparison of our approximate and exact results is given 
in sec. 3 and brief concluding remarks are presented in sec. 4. 



2 Theory 

2.1 Exact Lipkin-Model Calculation 

The Lipkin model consists of two single-particle levels labelled by a = — 
and +, each of which has a degeneracy p. We write the Hamiltonian 

H = Hq + V , where 
Ho = ^^J2 (^^la'^v'^ ' and 

pa 

pp'(j pp'cr 



2' 

pp'a 
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Here Hq is the unperturbed Hamiltonian with single particle energies ±2'^- 
The two-body interaction, V, has three terms. The interaction V acts be- 
tween a pair of particles with parallel spins and changes the spins from ++ to 

, or vice versa. The interaction is a spin-exchange interaction and U, 

which was not present in the original model flips the spin of one particle. 
It is of interest to note that the interaction does not change the value of the 
degeneracy labels pp'. 

Since each particle has only two possible states, the use of the quasi-spin 
formulation was suggested by Lipkin et al. [|l|. The quasi-spin operators 
obey angular momentum commutation relations and are defined by 

Jz = l^crCpg-apo- , J+ = gp+gp- , J_ = ^gp_gp+ . (2) 

pa p p 

The Hamiltonian can then be compactly expressed in the form 

H = eJ.+|X^(J+ + J') + |Vt^(J+J- + J-J+-n) + if/(J+ + J_)(n-l) , (3) 

where the number operator n = J2pa (^la^^p^ ■ The operator ,P = \{Jj^J^ + 
J-J+) + J1 commutes with the Hamiltonian so the Hamiltonian matrix breaks 
up into submatrices of dimension 2J+1, each associated with different values 
of J; for a given number of particles the largest angular momentum corre- 
sponds to J = |A^. It is straightforward to use standard angular momentum 
techniques to set up these submatrices which can then be diagonalized. Using 
a label a to distinguish the eigenvalues e, we have e = e(A^, J, a). 

Having obtained the exact eigenvalues for all A^, we can calculate the 
grand potential corresponding to the grand canonical ensemble according to 
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Q = —(3 ^InZ , where 

Z=J2 d{N, j)e-/^{-(^.^.")-/^^} . (4) 

NJa 

Here /i is the chemical potential and P = is the inverse temperature. 
The quantity d{N,J) gives the degeneracy, i.e., the number of times the 
angular momentum J occurs for a given A^-particle system. The physical 
quantities of interest, namely the mean number of particles (A^) and the 
internal energy E, can then be obtained from the grand potential with the 
usual thermodynamic relations 
dQ d 

(N) = -^ ' E=^m+^^{N). (5) 
This yields 

NJa 

E = Z-^Y1 e(A^, oi)d{N, j)e-/3{-(^.^'")-/^^} . (6) 

NJa 

2.2 Many-Body Approximations 

Our basic many-body approach is the Hartree-Fock approximation for which 
the finite temperature formalism is well known j^j. The HF single particle 
equations are 

eAk = {k\H,\i) + Y.{kj\V\i3)f, , (7) 
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where HF eigenvectors are used, i.e., \i) = q(+)|+) + Cj(— )|— ). Also the 
Fermi occupation probabihties are fj = [1 + exp(/9ej)]^^ with the definition 
ij = ej — //. The grand potential is then 

^^HF = -r' (l + e-^'') - E(u|V|zj)/./, , (8) 

i i>j 

from which, using eq. (5), the standard relations follow 

i 

E^^ = Y^e,U-Y.{^j\V\ij)M,. (9) 

i i>j 

Now we also want the grand potential giving the sum of the particle- 
particle hole-hole ring diagrams. This was evaluated in refs. 0, |[, but 
actually a much simpler derivation can be given in just a few lines. It is 
worthwhile to present this here. The ring series takes the form |^, [l^ 

Oring = f3-' E e^'-^^^TriFV - ^{FVf + liFVf - . . .] 

= (3-^ E e'^^'^^Tr ln(l + FV) , (10) 

where the notation FV means Fjj(zco'j^)(zj|V|A;/) and the summation is over 
i > j, etc. Here the Matsubara frequency = 2niyT, with u running over 
all integers. It should also be understood that a HF basis is used to evaluate 
the various quantities which arise. The pair propagator 

F.A^u^) = — , (11) 
lUi, - eij 

where we have introduced the simplifying notation Q{ij) = (1 — /i)(l — fj) — 
fifj = 1 — fi — fj and iij = ii + ij. Using a matrix notation, eq. (10) can be 
written 



The sum over v can be performed using eq. (13) of ref. |^ with the result 

-g-/3(el+QV) _ I 



^ring = /3 ""^Trln 



e-/3fi _ 1 

Choosing a diagonal representation, viz. 



(13) 



Y.{^^AoM + Q{i3){i3\V\kl)}m^n) = An{^3\Xn) , (14) 

k>l 

and noting that Tr In = In det, we obtain the final form 

Thus the grand potential is the difference between that obtained with ran- 
dom phase approximation (RPA) bosons and that obtained with unperturbed 
fermion pairs, treated as bosons. Now we should not simply add the ring and 
the HF results because the first order term of eq. (10) has already been in- 
cluded in the HF contribution. Thus we must subtract this, taking 

^^nng = ^^nng-E(^^1^K^')/^/. • (16) 

i>j 

Then fitotal = ^HF + ^ring- 

A simple, but approximate, method of obtaining the thermodynamic 
quantities of interest is to use the HF result (9) for the average number 
of particles and take the total energy to be -Etotai = -Erf + ^ring- This ignores 
the effect of the derivatives of eq. (5) upon fi^j^g and is therefore easy to com- 
pute. We shall comment upon this approximation later. Our aim, however, 
is to evaluate (A^)totai and i?totai exactly within the HF RPA formalism and 
this requires the derivatives of the HF and RPA energies and the derivative 
of the HF wavefunctions. 
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2.2.1 Evaluation of Derivatives 

Consider a general eigenvalue equation for a matrix which may, in general, 
be non-symmetric 

Cj\H\i)^EAj with {3\i)^5ij, (17) 

where the vectors \i) are the biorthogonal complements to the vectors 
Denoting partial derivatives with respect to some thermodynamic variable, 
X, by a prime, we have 

{y\H\i) + Q\H'\{) + Q\H\i) = E% , (18) 

6"K) + QV) = . (19) 

In the diagonal case, j = i, these equations give 

{i\H'\i) = El , (20) 
and in the off-diagonal case, j ^ i, 

{jm)^{E,-E,)m. (21) 

For the HF case we are dealing with a symmetric matrix so that \i) = 
Then applying the above to eq. (7), we have 

e'iS^k = (6.-6,)(^K')+E [{kj'\m)fj + {kjlVmfj + {kmtj)f'j\ .(22) 

j 

The derivative /j = ^ = -fj{l - fj)^\l3{ej - /x)], which involves the 
unknown e'y The derivatives of the HF wavefunctions are constrained by eq. 
(19) which implies — and — In our case there are only 

two states and labeUing these |1) and |2), we have 

|l') = c^|2) , |2') = -rf|l), (23) 



where d is a constant. Thus we have three unknowns €[, e'2 and d and these 
can be obtained by solving the three independant equations (22). It is also 
useful to note that eq. (22) yields the relation 



i>j 



(24) 



We also need the derivatives of the RPA eigenvalues A„. Using the vectors 
X which are biorthogonal to the vectors X, namely 



^{Xrn\ij){ij\Xn) = Smn , 
i>j 

it follows from eqs. (14) and (20) that we can write 



dx 



i>j,k>i ^ ^-^ 



+ Q{ij){ij\V\kl) 



(25) 



{kl\Xn) . (26) 



The derivative of the quantity in square brackets involves e^, // and which, 
as we have discussed, are obtained from eq. (22). 

Using eq. (5) and the grand potential Jltotai, we then have the exact 
expression for the number of particles 

-1 dA„ 



(TV) 



total 



I) -Eh--! 



i>j 



dfj, 

-1 diij^ 
dfj, 



(27) 



The internal energy is given by 



OK 
dp 



i>j 



d~€ 



dp 



-2J2{^m^J)f^fJ, (28) 



where I ^ - 



3 Calculation and Results 



The first step in carrying out the calculations is to compute the two HF 
energies via eq. (7) which requires that both the wave function amplitudes 
Cj(±) and the occupation probabilities fi be self-consistent. The derivatives 
of the HF eigenvalues and eigenfunctions can then be obtained from eq. (22) 
and the RPA eigenvalue equation (14) solved. Now the chemical potential fj, 
must be chosen to reproduce the correct number of particles. Initially this 
is done in the HF loop using eq. (9), but once faring has been computed the 
number of particles must be obtained from eq. (27). This will no longer be 
the desired value, so /i has to be adjusted and another iteration carried out 
and so on until the correct number of particles is obtained. 

For the RPA equation (14) we need to consider a five-dimensional basis, 
since the results are independant of the degeneracy labels p,p'. Specifically 
the basis is 

42^1 (29) 

where 1 and 2 label the HF states. The RPA matrix actually breaks into 
a 4 X 4 matrix (for p ^ p') and a 1 x 1 matrix (p = p') . Although we did 
not use it to simplify the calculations, we point out that in actuality the 
4x4 matrix can be split into a 1 x 1 matrix corresponding to the linear 
combination {o^pi20^pi ~ o^pii0^p2j ^"^^ a 3 x 3 matrix for the orthogonal states. 

We shall discuss the case where the available states are half-filled, i.e., 
(N) — p, since other choices do not yield qualitatively different results. Also 
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for the pure HF case the levels 1 and 2 must either be completely filled or 
completely empty at T = 0, so the choice (N) = p yields a well-defined 
T = limit. For the half- filled case the pure HF result from eq. (9) requires 
a chemical potential // = |(ei + £2), which implies that /i + /2 = 1 or 
Q(12) = (5(21) = 0. This means that only the a^i-^^ali and 0^,20^2 states yield 
a non-zero RPA contribution. When the chemical potential is obtained from 
eq. (27) which employs fitotai, this is no longer precisely true, neverthelesss 
the contribution of the three "12" states remains small. 

We shall present results for the case where p = 16 and (A^) = 16 so that 
half of the 32 available states are filled. We have examined other values of 
(N) = p and found no qualitative differences, although the approximations 
are quantitatively a little less accurate for smaller numbers of particles, as 
one might expect. We will take 1^ = 1, thus implicitly measuring energies in 
units of ^, i.e., the quantities we discuss are dimensionless. Further we will 
choose U = W since no qualitative difference is observed if they are unequal. 

For our first set of calculations we chose for the parameters of the Hamil- 
tonian of eq. (1) U = W = and V = —0.65 so that we can discuss the 
situation addressed in refs. 0, ^, ^ . With only V non-zero analytical solu- 
tions can be obtained for the pure HF case 0, 0. Defining ei = —62 = — |e^^, 
a solution with e^^ = 1 is always possible, i.e., in this case the HF Hamil- 
tonian is unchanged from the unperturbed Hamiltonian, Hq. However, if it 
exists, the solution of 

e^F^ |y|^p_l)tanhi/3e"^ (30) 

gives a lower energy [@, ^. As T increases, i.e., (3 decreases, the tanh de- 
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creases until the limiting case is reached where e = 1; for our parameter 
choice the critical temperature = 2.43. This behavior is illustrated by the 
dashed curve in fig. 1 where we plot the lower eigenvalue ei. The numbers on 
the curves indicate the intensity Ci(— )^. This starts out at 0.55 and becomes 
unity after the phase transition to the unperturbed state. The corresponding 
internal energy, -Erf, is shown as a function of temperature in fig. 2 (dashed 
curve). At the phase transition the slope becomes discontinuous and for 
higher temperatures the energy arises from the unperturbed Hamiltonian, 
Hq, only. This can be contrasted with the exact result denoted by diamonds 
in fig. 2 where the curve is smooth and there is no indication of a phase tran- 
sition. We conclude that the discontinuity is an artifact of the approximation 
which is employed, without physical significance. The effect of including the 
ring diagrams is indicated by the solid curve in fig. 2. We see that that they 
yield a modest improvement in the results at low temperatures, in agreement 
with the zero-temperature results of ref . . At fairly high temperatures the 
effect is larger and brings the calculations close to the exact result. However, 
the rings yield a discontinuous curve in the "phase transition" region and 
the results are inaccurate there. We remark in passing that a similar "phase 
transition" appears for the case U = V = with W ^ and again the exact 
calculations show no evidence for such an effect. 

We may note that, as T — > oo, the occupation probabilities /« — > | and 
therefore Q{ij) ^ so that the RPA energies A„ of eq. (14) are just the HF 
energies iij. It is straightforward to check that in this limit ^I'^^^g —>■ 0. The 
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internal energy is just the HF energy which can be written 

i?total = i?HF = i JZi^mtj) = -\PW . (31) 

Since the pair propagator vanishes {Q{ij) — > 0) and the particle-hole prop- 
agator also vanishes (because it is proportional to [(1 — fi)fj — fi{l — fj)] 
which is zero in this limit) this should be an exact result. That is, at infinite 
temperature only contributions of first order in V survive. We have verified 
that eq. (31) agrees with the results of our exact calculations. For the case 
shown in fig. 2, eq. (31) implies that the asymptotic internal energy is zero. 

We next examine the effect of taking the relatively modest values U = 
W — —0.02, with the same value of 1^ = —0.65. The results for this case are 
given in fig. 3. As compared with fig. 2, the change in the internal energy is 
small and the exact results are quite similar in the two cases. However there 
is a qualitative difference for the HF curve which now smoothly approaches 
the unperturbed result for high temperatures. The unperturbed case, which 
arises from Hq and is the same for all the calculations we present, is repre- 
sented by the dot-dashed curve in fig. 3. Comparing this to the other curves 
at low temperatures, we see that the effect of the perturbation V is very 
large indeed . The corresponding HF single particle energy and intensity in 
fig. 1 arc similar to before except, that there is no phase transition and that 
asymptotically for T ^ oo the intensity Ci(— )^ 0.98 rather than unity. 
The effect of including the ring diagrams here is shown by the sohd curve of 
fig. 3. This differs from the previous case in the region where the transition 
from the low to the high temperature behavior takes place- the curve is now 
smooth and the agreement with the exact results is much better. 
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As a final case we take U = W = —0.2, these values being comparable to 
V, which again is —0.65. The results in fig. 4 show that the internal energy, 
E, is roughly doubled at low temperature and the pure HF approximation 
gives very good agreement with the exact answer. The effect of the ring 
digrams is small, but they do provide even better results at high tempera- 
tures. Reducing the rather large value of V that we have used to make it 
comparable to or less than U and W results in even less of an effect from 
the ring diagrams. This is not unexpected since the matrix element between 
the 11 and 22 states dominates the RPA correlations and this is strongly 
influenced by V. The single particle energies here (solid curve of flg. 1) are 
much larger than in the previous examples and the mixing between the basis 
states remains large at high temperature. In fact for T ^ oo the intensity 
ci(-)2 ^ 0.66. 

Finally let us discuss the approximation of taking E = Euf + ^ring- 
T — > the quantity /j(l — fj) goes exponentially to zero, in which case 
eq. (22) indicates that the derivatives of the HF energies and wavefunctions 
become zero. Thus setting the derivative contributions in eq. (28) to zero, 
the T = ring contribution to the energy is just —(J2n — J2i>j ^ij), where 
the summation runs over those states for which A„ and iij < 0. In other 
words the states with energies less (greater) than the chemical potential are 
fllled (unflUed). The same expression is obtained directly from firing in eq. 
(15) and this is the well-known T = result 0. Thus, at T = 0, the 
internal energy is exactly given by E'hf + ^ring- However as the temperature 
increases this becomes an approximation and it begins to deteriorate when 
the slope of the internal energy curves in flgs. 2-4 starts to increase. Indeed 
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filing is always negative, whereas (E'totai — -E-hf) is positive in some regions. 
In the high temperature regime £^hf + ^Hng gives energies that are roughly 
halfway between the HF and the HF+ring results. Thus the accuracy of this 
approximation can only be relied upon at low temperatures. 

4 Concluding Remarks 

We have calculated the thermodynamic properties of the Lipkin model Hamil- 
tonian exactly and compared with approximate many-body treatments. The 
case where only the V (or only the W) term of the Hamiltonian is non-zero is 
special because the Hartrcc-Fock single particle states differ from the unper- 
turbed values only for temperatures up to some critical Tc. The transition 
between the two situations manifests itself as a discontinuity in the slope of 
the calculated quantities. However no such effect is observed in the exact 
calculations. Further if the other parameters of the Hamiltonian are allowed 
to differ from zero, even by a relatively small amount, this behavior of the 
HF theory disappears and the calculated curves are smooth with no disconti- 
nuities in the derivatives. For both of these reasons we conclude that the HF 
"phase transition" is an artifact which is not likely to be relevant to actual 
physics. 

In cases where all the parameters of the Hamiltonian arc non-zero, which 
we think is much more likely to be representative of actual situations, we 
find very good agreement between the exact and the HF results. We found 
this somewhat surprising since the differences between the exact and the 
unperturbed internal energies can be large, i.e., the perturbation is not in any 
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sense small. We also examined the effects of the particle-particle ring series 
and found that accuracy demands a thermodynamically correct treatment 
of this contribution, i.e., the number of particles and the internal energy 
should be obtained from the thermodynamic relations (5) using the complete 
grand potential Ototai- The ring effects yield a modest improvement of the 
HF result which is mainly evident at high temperatures, although there is 
a small effect in the low-T regime. Regarding the ring contribution, we 
make two cautionary remarks. Firstly we have used a rather large value of 
the parameter V so as to investigate the situation discussed in the previous 
paragraph; this has the effect of enhancing the size of the rings. Secondly in 
the infinite temperature limit we have pointed out that the ring contribution 
becomes zero. Nevertheless our full approximation is remarkably accurate 
and it would be interesting to see what results it yields for more realistic 
Hamiltonians. 
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Figure Captions 



Figure 1. The Hartree-Fock energy of the lower state, function of 

temperature for the three different parameter sets indicated. The numbers 
indicate the intensity of the |— ) component of the corresponding eigenvector, 
i.e., ci(-)2. 

Figure 2. Comparison of the HF, the HF and ring diagram and the exact 
values of the internal energy as a function of temperature. The three cases 
are denoted respectively by the dashed curve, solid curve and by diamonds. 
The parameters are V — —0.65, U = 1^ = 0. 

Figure 3. As for fig. 2, but with parameters V = -0.65, U ^ -0.02. 
Here the dot-dashed curve gives the result obtained with the unperturbed 

Hamiltonian, Hq, only. 

Figure 4. As for fig. 2, but with parameters V — —0.65, U — W — —0.2. 
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